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The structure of the jet in Cyg X-1 inferred from orbital modulation 
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ABSTRACT 

We study free-free absorption of radio emission by winds of massive stars. We derive formulae 
for the optical depth through the wind measured from a point of emission along a jet, taking 
into account Compton and photoionization heating and Compton, recombination, line and 
advection cooling. 

We apply the developed formalism to radio monitoring data for Cyg X-1, which allows 
us to obtain strong constraints on the structure of its inner jet. With the data at 15 GHz, and 
taking into account an anisotropy of the stellar wind in Cyg X-1, we estimate the location of 
the peak of that emission along the jet at about one orbital separation, i.e., ~ 3x10'^ cm. Given 
a previous determination of the turnover frequency in Cyg X- 1 , this implies the location of the 
base of the jet at ~ 10^ gravitational radii. We also obtain corresponding results at 8.3 GHz 
and 2.25 GHz, which roughly follow the standard conical partially self-absorbed jet model. 
Furthermore, we find that the level of the orbital modulation depends on the radio flux, with 
the modulation being substantially stronger when the radio flux is lower This is explained 
by the height of the radio emission along the jet decreasing with the decreasing radio flux, as 
predicted by jet models. Based on the finding of the flux-dependent orbital modulation, we 
are able to estimate a range of the possible changes of the form of the radio/X-ray correlation 
in Cyg X-1 due to free-free absorption. We also derive predictions for the orbital modulation 
and flux attenuation at frequencies beyond the 2.25-15 GHz range. 

Key words: accretion, accretion discs - radio continuum: stars - stars: individual: Cyg X-1 
- stars: individual: HDE 226868 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 



Cyg X-1 is an archetypical and well studied persistent black-hole 
binary. Its companion is the OB supergiant HDE 226868. Cyg 
X-1 is both a radio and X-ray source, and the fluxes in these 
two bands are strongly corr elated (e.g., lOallo. Fender & PoolevI 
l2003l : IZdziarskietani2011bl. hereafter Paper I). The likely cause 
of the correlation is the jet in this system (Stirling et al. 200 ll . here- 
after SOI) being formed by the matter of an inner hot accretion 
flow iHeinz & Sunvaevl l2003l : iMerloni. Heinz & di Matted [20031 : 
lYuan & Cui 2005^ The radio emission is modulated at the orbital 
period (Poole v. Fender & Brocksopnl Il999l) , which is caused by 
free-free absorpti on in the stellar wind from the massive compan- 
ion. As shown bv lSzostek & Zdziarskil ( [2007h (hereafter SZ07), the 
observed modulation levels imply that the average radio fluxes are 
substantially absorbed, which then, as discussed in Paper I, modi- 
fies the form of the observed radio/X-ray correlation. 

In this work, we study in detail the free-free absorption of ra- 
dio emission in winds of massive stars in high-mass X-ray binaries 
(hereafter HMXBs). We derive formulae for the optical depth mea- 



sured from the point of emission along the jet. We also consider 
emission of counter-jets and Compton scattering by the wind. 

Then, we apply our results to Cyg X- 1 , t aking into account the 
strong anisotropy of its stellar wind, found bv lOies et al.lfcOOSh . as 
well as the radio emission being extended (SOI). This yields strong 
constraints on the structure of the jet in Cyg X-1, in particular it 
implies that a major part of this emission takes place at distances 
comparable to the binary orbit. Based on the results of this applica- 
tion, we estimate the effect of the free-free absorption on the form 
of the radio/X-ray correlation in Cyg X-1. 



2 RADIATIVE HEATING AND FREE-FREE 
ABSORPTION IN A STELLAR WIND 

A standard wind velocity profile from a massive star is (e.g., 
iLamers. Cerruti-Sola & Perinotta,1987i) . 



v(r) ^ 



(1) 



where r is the distance from the centre of the donor of radius R,, Voo 
is the terminal velocity, and p parameterizes the wind acceleration. 
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Here, a small correction due to velocity reaching the sound speed 
rather than being null at the stellar surface has been neglected. The 
ion density is. 



n(r) 



(2) 



where / is the direction-dependent factor accounting for the wind 
asymmetry due to its focusing (see below), M is the total mass- 
loss rate, & and <p (measured from the axis connecting the stars) are 
the polar and azimuthal angle, respectively, is the proton mass, 
/Ji - 4/(1 + 3X) is the mean ion molecular weight, and X is the H 
fraction. For simplicity, we hereafter assume X = I. 

The local wind temperature follows from the balance of radia- 
tive heat ing and cooling and advecti on. In an optically-thin approx- 
imation jKallman & McCravlll98 j) . applicable here, the radiative 
rates depend only on the shape of the irradiating spectrum and on 
the ionization parameter, 

f - (3) 



dTarter. Tucker & SalDetejri969h . where x is the distance to the ir- 
radiating source with the ionizing luminosity, Li„„. The net rate of 
X-ray photoioniz ation heating a nd recombination cooling per ion 
has been fitted as telondiiJll994h . 



Cx ^ 1.5 X 10"^'^"'*r-"'(l - T/Tx)n 



(4) 



(in units of erg s"'), where T is in K and kTx - 10 keV is the 
characteristic energy of the ionizing spectrum. Given that the wind 
temperature in Cyg X-1 has kT ^ I keV, we can neglect t he sec- 
ond term in parentheses above. The formulae of Blondinl ( Il994h 
have been derived for the case of a 10-keV bremsstrahlung X-ray 
spectrum in a high-mass X-ray bi nary. The rate o f bremsstrahlung 
and line cooling has been fitted as ( lBlondiiJll994h . 



-Cb,i ^ 3.3x 10-"'r"Xr)+ 

[l.7 X 10^'* exp(-TJT)^-'T-^'^ + 10"^"] n(r) 



(5) 



where = 1.3 x 10^ K. 

In the non-relativistic regime, the rate of Compton energy ex- 
change between electrons and photons is proportional to j(E - 
4kT)F(E)dE. The rate per electron including Klein-Nishina cor- 
rections but assuming kT <k m^c- can be written as, 

4kTa-T 



Cc = 



8(6)EF(E)dE ■ 



h(<E)F{E)dE, 



(6) 



where F(E) is the differential energy flux, ctj is the Thomson 
cross section, is the electron mass, e = Elm^c^, and g{e) 
and h{e) give the Klein-Nishina corrections to C ompton heating 
and cooling, respectively. The f ormer is given by (lGuilbertl ll986l- 
[Beloborodov & Illarionovll T995l) . 



g{e) = 



3(6- - 26 - 



3)log(l +26) _^ 9 + 5Ie + 93e2 +51e^ - 106"* 



I 



2l6 



86^ 

1476^ 



46^(1 + 26)3 



+ 0(6^). 
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Separating F(E) into the stellar and X-ray component, we have. 



(7) 



g{e)ELx{E)dE + 



TjkT.L, - kT(L, + Lx) 



,(8) 



where L, is the stellar luminosity, T, is the stellar temperature, 
Lx(E) is the differential luminosity originating at the black-hole 
location, r and x is the distance to the centre of the donor and to the 
compact object, respectively, 77 = 90^(5)/;r"* ^ 0.96, and f is the 




Figure 1. The function g(E/m,:C^), by which the Compton heating of elec- 
trons with kT <s njef- by photons with energy E is reduced due to the 
Klein-Nishina effects. 



Riemann zeta function. The X-ray heating term in equation (|8) is 
valid at any values of 6, provided kT jm^c^ <k 1. Equation l[8j also 
assumes that cooling is dominated by photons at e « 1, which, 
e.g., is the case when L, » Lx- I f it is not satisfied, h{e) needs to b e 
included as well, see eq. (17) in lBeloborodov & Illarionovl l l 19951) . 
The function g{e) is plotted in Fig.[T] We see it decreases rather fast 
at £ > 50 keV. At a point (r, j?, ip). 



X = r + a~ - 2ar sin i? cos ip. 



(9) 



where a is the separation between the stars. 

Then, assuming the velocity profile of equation ([TJ, we need to 
solve the energy equation, which, for an ionized gas, can be written 
as. 



v(r) 



dr d 1 

2,k—+2kTn{r) — — 
dr dr n{r) 



■ Cx + Cb.i + Cc 



(10) 



This can be solved for T{r) adopting a boundary condition close 
to the stellar surface, where the wind is optically thick and our 
assu mption of optical thinness bre aks down. In an optically thick 
wind. lCarciofi & Biorkma"nl(l200d) find T ~ 0.6T,. We have found 
that the solution above the assumed boundary, r > ro, is virtually in- 
dependent of the value of ro, and we adopt ro = (1.2-1.3)^*. At low 
radii, dT/dr > 0, and cooler particles are brought to a given point 
from radii closer to the star At all radii, adiabatic expansion cools 
the wind. Thus, the dominant effect of the left-hand side of equation 
jlOl is cooling. We have found that solutions of this equation are 
rather well approximated by replacing the differential terms (i.e., its 
left- hand side) by an algebraic adv e ction cooling term, -Cadv, see, 
e.g., IChen. Abramowicz & Lasotal ( Il997h for a similar approach. 
We have found that we achieve the highest accuracy in this ap- 
proach defining C^dv as. 



3kTv„ 



We thus solve an algebraic equation for T(r), 



Cx + Cb.l + Cc + Ca, 



: 0. 



(11) 



(12) 



See Section|5]for an example comparison of the two solutions. 

We then assume the bulk of the radio emission at a given 
V occurs at a distance, z, along a jet from its origin at the 
compact object, see Fig. |2] This z corresponds to the place 
along the jet where it becomes optically thin to synchrotron self- 
absorption CBlandford & K onigl 1979, hereafter BK79). The jet 
is assumed to be perpendicular to the circular orbital plane, see 
Fig. |2] We note that the assumption of perpendicularity is not 
a crucial; the jet can be inclined (as suggested for Cyg X-I by 
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Figure 2. A schematic representation of the considered system at the su- 
perior conjunction. The distance between the high-mass donor (wliich cen- 
tre is at 0) and the compact object is a. The point of radio emission is at 
the height z along the jet (shown by the red hne), which originates at the 
compact-object location. A point along the photon path through the wind 
towards the observed is at the distance, /, from the emission point, at the 
distance, r, from the centre of the donor, and at the distance, x, from the 
compact object, which is surrounded by the source of X-rays. 



iMalzac. Belmont & Fabiajl2009h or culled (as in the Cyg X-1 
model of SZ07), which introduces a phase lag to an orbital modu- 
lation. Then, z/a is to be interpreted as the height above the orbital 
plane in units of the distance to the rotation axis of the donor. In 
this case, a would be somewhat higher than the orbital separation. 
We neglect here this complication. We also neglect any effects of 
the jet opening angle, as a second-order effect. In the case of Cyg 
X-1, it has been constrained to be £ 2° (SOI), which is negligibly 
small for our calculations. 

As shown in Fig. |2l / is the path length from the emission 
location along the jet to a point at a radius, r, from the donor centre. 
We find, 

(r/a)- = 1 + izlaf + (l/af - 2(1 /a) [sin ; cos (/> - (z/a) cos /] , (1 3) 

where / is the binary inclination and ip is the orbital phase, which 
equals tp = 0, n al the superior and inferior conjunction, respec- 
tively. Then, 

(x/af = (z/af + (I /a)- + 2(z/a)(l/a) cos i. (14) 
The free-free absorption coefficient at a frequency, y, is, 

Q-ff = O.llT-^'^n^v-- cm"' cm"', (15) 

where T is in units of K and v is in GHz. At a reference temperature, 
T = Tq, we can write ajf = Q'o(«i/«o)^. where Ug is the absorption 
coefficient at an ion density of iiq. We define ng as the density at r = 
a under the assumption of i' = Vc . We also define a normalization 
optical depth of Tq = aga. 

Then, the free-free optical depth is, 



To 



fire 4) 



T(r) 



To 



-3/2 



d(l/a). 



(16) 



This can be calculated numerically with T(r) as the solution of ei- 
ther equation dlOt or ( I12t . 



On the other hand, we also consider the case of an isothermal 
wind with constant velocity (which is approached at at r » R,). In 
that case. 



T - To 



d(l/a). 



(17) 



In this approximation, we obtain an analytical formula for t, 
T nt/2 + su^'^ -I- farctan(™""^) _ ?arccot(iM""^) - 

~, ~ 2?t(3/2 " 2?m3/2 ' 

t=\+(z/a)^, i = sin i cos - (z/a) cos I, u = t-s^. (18) 

We find « > everywhere except for ifi = n, z/a = cot /, when u = 
and equation l llSt appears to have a singularity. However, t is still 
finite there, t/tq = sin' i/3. For z —> 0, 



T n + 2iu + sin 2yu 



+ 0(z,i,(f>), ;u = arcsin(sin I cos 0), (19) 



To 4 cos^^ /J 

At t((I> = 0) - T(4> = 7T) 2[+sin2/ 

— = = T— + 0(z, i), (20) 

To To 2 COS^ I 

where At is the difference in the optical depths at the superior and 
inferior conjunctions. For large z, the leading order of t is indepen- 
dent of (f), 

T 2i - sin 2i . 

- = 3 + 0[(z/ar\ i, 0], z/a » 1, (21) 
To 4(z/a)3 sin ; 

and At is lower by one order of a/z, 

At 9 sin 12/ cos ! -I- sin 3; 

To 4(z/a)'* sin"* ! 



+ 0[(z/ar\i], z/a»l. (22) 



Potential constraints from eclipses of the radio emission from 
both the jet and counter-jet are discussed in Appendix lAl Appendix 
|B]presents results for Compton scattering (or another process with 
the cross section proportional to the density) analogous to the above 
ones for free-free absorption. 

We define the full modulation depth due to absorption as. 



„ ^(<^max) - /^(^ram) , , . , 

D = TT—. ^ = 1 - exp(-AT), 



(23) 



/^(f^max) 

where F((^n,ax) and F(0„,in) are the maximum and minimum fluxes 
in the model averaged and folded light curves, observed at (^„,i,x and 
(/>mm, respectively, and 0,ni„ represents the possible phase lag. In our 
model, (^max - 0min = ^- We note that D (or, equivalently. At) in 
the model is uniquely determined by z/a, i and To. The attenuation 
averaged over the orbit and the average t are defined as, 

r exp(-T)d0 

A = = exp(-f), (24) 

where the integration can be done over only 0-n due to the planar 
symmetry. As a rough approximation for t, T(n/2) may be used. 
We note that D, A, f , 0n,i„ are functions of frequency. 



3 PARAMETERS OF THE CYG X-1 BINARY AND ITS 
STELLAR WIND 

The ephemeris of Cyg X-1 has been determined, e.g., by 
iLaSala et a i]( ll998h . lBro'cksopp et alJ ( fT99^ . It can be written as. 



„[MJD] = 50077.995 + PE, 



(25) 



where P = 5.599829 d, t^up is the time of a superior conjunction 
(the black hole furthest from the observer), and E is an integer. 

Observational determinations of the mass of the donor, M,, 
and of the black hole, Mx, R,, i, and the distance, d, are all mutually 
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conne cted. We use here the constraints of 'Caballe ro-Nieves et alj 
1 20091). and of the evolution ary models of iZiolkowski 1 2005h . 



ICaballero-Nieves et al.l ( l2009t) obtained the most likely ranges of 
M, ^ (17-31)Mo, Mx ^ (8-16)Mo, R, ^ (15-21)%, / ^ 43°- 
31°, d - (1.5-2.0) kpc. They also measured the donor temperature 
as T, = 2.8 X lO"* K, and derived. 



RJRq = (10.3±0.7)d, 



:arcsin[1.02±0.09)/rf], 



(26) 



where d is in kpc. Here, we choose their sol ution that also agree s 
with the detailed evolutionary modelling of IZiolkowski I ( |2005|) . 
namely, ( = 31°, d = 2.0 kpc, M, = 30Mo, Mx = 16Mq, 
R, = 20.7Ro, which yields a 47.5Ro - 2.3i?J3 The eccen- 
tricity is <K 1 ( Orosz et al. 201 1), and we here assume the orbit to 
be circular. 

The total ma ss loss rate in the hard spectral state was found by 
iGies et al.1 i2003l) as -M ^ 2.6 X IO^'^Mq yr"'. However, the mass 
loss can be highly asymmetric in binary systems, especially those in 
which the donor is close to filling its Roche lobe. T his is the case in 
Cyg X -1 , where R, > 0.9 of the Roche-lobe radius dOies & BoltonI 
Il986al lbl: loies et alj[2003l : IZiolkowski 1120051). In this situation, th e 
wind is focused towards the black hole Ipriend & Castoijll982l) . 
with the wind density inside the Roche lobe being substantially 
higher than that outside it. The presence of the focused wind com- 
ponent was found to fit very well the He n /14686 emission profile 
of Cyg X-1 is les & Boltonil986hl) . Evidence for the focused com- 
ponent bas ed on soft X-ray li nes was also found by iMiller et al.l 
(12005) and Ha nke etal] j2009t) . 

Further com plexity of the wind geometry was found by 
lOies et all j2008l) . Namely, they found the unfocused, locally spher- 
ically symmetric wind component emitted by the X-ray irradiated 
hemisphere of the donor to be undetectable. This weakness is ap- 
parently due to wind ionization by X-rays, which strongly reduces 
the wind radiative acceleration by the lines emitted by the donor. 
On the other hand, the other hemisphere of the donor, shaded from 
the X-rays, emits a normal wind, similar to that of an isolated mas- 
sive star (Blondin 1994). Summarizing, the wind of HDE 226868 
has three main components: the wind focused along the axis joining 
the stars, which carries most of the mass loss; the standard radia- 
tively driven wind from the hemisphere in the shadow, and a weak 
wind component from the donor hemisphere facing the black hole. 
As found by Gies elak ( 2008), this third component is weak in both 
spectral states, hard and soft, of Cyg X-1. 

Radio photons are emitted by the jet relatively high above the 
orbital plane and may be attenuated mainly by this third weak com- 
ponent, i.e., the quasi-spherical wind from the hemisphere facing 
the X-ray source. Given its weakness, it has / <K 1, but no specific 
constraints are currently available. For the wind velocity profile, we 
use here t hat of equation ([TJ with i'„ 1.6 x 10^ cm s"' and jS = I 
dOies & B olton 1986b). 

The stellar luminosity for our assumed parameters is L, ^ 
9.1 X 10^^** erg s"'. We then use the average hard X-ray/soft 
y-ray spectrum of Cyg X-1 as observed by the OSSE and 
BATSE detectors on board of Compton Gamma Ray Observatory 



' After the calculations reported where completed, iReid et al] \2Q\ lb have 
determined the trigonometric parallax distance as d = 1.86|'^jJ"[j kpc (with 
Itr unce rtainties), wh i ch val ue is compatible with adopted here within ~ 
\<T. Then. lbrosz et aljbOl ll) determined My = 14.8 ±1.0Mo, again within 
Icr from our value. They also found M, = 19.2 ± 1.9M0, R, = 16.2 ± 
0.7Rq, and i = 27.1 ±0.8°. The effect of these difference on our results have 
been found to be relatively minor, and less important than the uncertainty 
regarding the structure of the stellar wind. 



dMcConnell et al.1 l2002h . We fit this spectrum with an e-folded 
power law, Fx(E) = K(E/lkeYy exp(-E/Ef ), obtaining K = 
0.62 cm S-' (normalized to the OSSE flux), a = 0.36 ± 0.03, 
Ef = 160 ± 5 keV, yielding Lx ^ 1-7 x 10" erg s"'. The fit is 
based on the data at > 20 keV, and it quite severely underestimates 
the actual spectrum in softer X-rays. However, the main contribu- 
tion to Compton heating is from hard X-rays, see equation (|8), 
and thus this underestimate has a negligible effect. We note that 
the Klein-Nishina correction is very important for this spectrum, 
with the Compton heating rate calculated non-relativistically, with 
g(6) = 1, being 2.7 times higher than the correct rate calculated 
using g(6) of equation ([7}. On the other hand, photoionization is 
mostly due to soft X-rays, which spectrum in Cyg X-1 is relatively 
uncertain due to absorption. We assume here Lj^n = 1 x 10^^' erg s"' 
in equation 

Then, we can express the characteristic free-free optical depth 
of the wind, To in equation J16l >. as. 



To ^ 35.5 



-fM 



2.6 X lO-'^M© yr- 



6 X 10* cm s" 



\10«K/ 
1/ V 15 GHz/ 



-3/2 



M, + Mx] 
46Mo j 



(27) 



where we used 10* K as the reference temperature. Giv en the lack 
of a measurement of the wind in the polar regions dCies et al.l 
|2008|) . we treat / as a free parameter of the model. 



4 RADIO EMISSION IN THE HARD STATE OF CYG X-I 

We study here radio monitoring data at 15 GHz from the 
Ryle Telescope, which cover MJD 50226-53902, and the Ar- 
cminute Microkelvin Imager (AMI), MJD 54573-55540. Also, 
we use the 2.25 GHz and 8.3 GHz monitoring data from 
the Green Bank Interferometer, for MJD 50409-51823. These 
data are the same as in Paper I, and we refer to it for de- 
tails. We fit only data corresponding to the hard spectral state, 
which corresponds to MJD 50350-50590, 50660-50995, 51025- 
51400, 51640-51840, 51960-52100, 52565-52770, 52880-52975, 
53115-53174, 53540-53800, 53880-55375 (which intervals con- 
tain 12408 Ryle/AMI measurements a nd 1160 GBI ones), see 
IPoutanen. Zdziarski &Ibragimovld2008b and Paper I. 

An important issue here is that these telescopes do not spa- 
tially resolve the radio emission. On the other hand, the emission in 
the hard state is known to be extended. SOI have found that ~ 1/3- 
1/2 of emission at 8.4 GHz observed by VLBA is within the beam 
of a longitudinal extent of ~ 3", whereas the remainder forms an 
extended jet visible up to ~ 15". This has been confirmed by further 
VLBA observations at 8.4 GHz reported in Rushton et al. (201ll), 
in which the resolved fraction was ~ 0.3-0.4. The resolution of 3" 
corresponds to the length of 9 x 10'''(d/2, kpc)/ sin / cm, which is 
~ 30fl/ sin i. The resolved 8 GHz emission can thus have only a 
tiny orbital modulation. 

The observational situation at 15 GHz is mo re complex. The 
VLBA beam size of the observations reported in iRushtonI d2009h 
is ~ 1.3", which corresponds to ^ 4 x 10'^((i/2, kpc)/ sini cm 
10a/ sin The unresolved and total fluxes for the three observa- 
tions taken on the same days as the 8.4 GHz ones of SOI (i.e., MJD 
5 1035, 5 1037, 5 1039) are (7.3, 7.7) mjy, (10.9, 12.2) mjy, (6.1, 7.8) 
mjy. Thus, the unresolved fractions are relatively small, 0.05-0.22. 
However, Cyg X- 1 was also monitored by the Ryle Telescope three 
times on each of those days, yielding the total fluxes of (13.2, 14.9, 
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16.0 mJy), (13.2, 12.2, 10.4 mjy), (18.9, 19.4, 18.1 mjy), respec- 
tively, i.e., substantially more than the total VLBA fluxes for the 
1st and 3rd observation. Although all of the Ryle observations were 
taken ~ 0.1 d before the start of the VLBA observations, Cyg X-1 
was in a relatively stable state at that time with no large-amplitude 
variability detected. It thus appears that the VLBA measurement 
missed an extended part of the resolved emission. Then, the re- 
solved fractions are ^ 0.1-0.7. 

As discussed in lZdziarski. Lubinski & Sikoral ( l2012l) . It IS un- 
likely that the radio jet in Cyg X-1 follows a single set of power-law 
dependencies from its base to the resolved parts. It appears that it 
contains a compact part within ~ (l-lO)a, and then a secondary 
dissipation event leads to formation of a more extended part of 
the jet, resolved by VLBA. The observed strong orbital modula- 
tion (see below) has to occur in the compact part. The division of 
the fluxes between the two parts as a function of time remains un- 
known. Based on the above data, we can only roughly estimate that 
about a half of the flux is emitted by the compact part. Given that 
the radio-IR spectrum of Cyg X- 1 appe ars to be a single power law 
jpender et al.l200(]| : lRahoui et alj20n"h . both components have ap- 
proximately the same spectra. 



5 FREE-FREE ABSORPTION IN CYG X-1 

With the model developed in Section|2l we can fit the orbital mod- 
ulation profiles, which are given by the light curves averaged and 
folded over the orbital period. In this way, we obtain the location 
of the bulk of the emission, z/a. We neglect the emission of the 
counter-jet, as its emission is beamed away from the observer and 
it passes through a much higher column density than that of the jet, 
including the dense focused component close to the binary plane. 
Given the results discussed in Section |4] we split the radio light 
curves into unmodulated and modulated parts. 



F{t) - fc-F(f)iinraodulatcd + (1 " i)^^(f)raodulatcd , 



(28) 



where t it time. The second term on the right-hand side corresponds 
to emission of the inner part of the jet, within ~ 10' a or so, which is 
orbitally modulated. The first term corresponds to the larger-scale 
emission, which is not modulated. Then, 



Fi„ir(0 = f(0[& + (l-*)exp(f)], 



(29) 



where F,ntr, F denote the light curves averaged over time intervals 
of the length P. Given that the obs erved emission comes from a 
regio n of the size < 10'^ cm (SOI: iRush ton 2009; Rushton et al. 
l201lh and a relativistic speed of the jet (SOI; .Malzac et al.,,2009,) . 
we assume that both parts of Fj,,,,. vary in the same way (apart of the 
orbital modulation) on a ^ lO'^'^-s time scale, below which scale the 
radio variability is weak. 

Equation ( 1291 ) implies the average attenuation of 



1 



(30) 



D = 



(31) 



b + {\- Z7)exp(-f) 

The modulation depth becomes, 

1 - exp(-AT) 
1 +Z?exp(r„,in -f)/(l -fo)' 

where T„,in is the minimum optical depth through the wind (corre- 
sponding to the maximum flux, at ^mzud- 

In the calculations below, we assume either b = 0.5 (see Sec- 
tion|4]l ovb = 0, which provide us with an estimate of the sensitivity 
of the calculations to taking into account the unmodulated part of 



the flux. The assumed b = 0.5 is, given the available data, only 
approximate, as well as it is unlikely that it was constant during 
the course of the monitoring. In Figs. |3}j6]below, we show the re- 
sults for b = 0, which directly correspond to the observations of 
the modulation. After obtaining folded and averaged light curves 
of F(t), we subtract from them b{F), as implied by equation l |28l l, 
and then recalculate the folded average modulation. 

We need to properly average the folded light curve. We are in- 
terested in variability governed by the orbital time scale, and not in 
that on much shorter time scales. Therefore, we rebin the observed 
light curve with a bin size equal AP = P/K, where K is the num- 
ber of bins per period. In this way, we also partly remove a bias 
on the folded averages due to non-unifoiTn coverage, e.g., due to 
a large numbers of points in some orbital bins and low number in 
other bins. The errors on the resulting averages are then calculated 
from dispersion of fluxes averaged within individual bins of the 
same phase range and the size AP. This method was applied to pe- 
riodic variability of Cyg X-1 by Ibragimov, Zdziarski & Poutane^ 
( l2007h . We fit J^^ exp(-T)d(f)/A(p, where Ac/> = In IK is the phase 
bin width, instead of using the value at a bin centre. Also, finite 
phase lags with respect to the superior conjunction are observed in 
Cyg X-1. Then, the phase lag is a free parameter of the model, and 
t(0) of equations l ll61418t is replaced by t(0 - ^min)- 

We first compare our results with those of SZ07 (whose re- 
sults are for b = 0). In that work, the atomic heating and cooling 
was treated more accurately than in this work. On the other hand, 
adiabatic cooling was not included, which resulted in a constant 
temperature, Tg, at large radii. Thus, we first compare their results 
with those using our analytical, constant T, model, equation il8t . 
Using the parameters from SZ07, namely, f = l,T = 3.2 x 10^ K 
and I = 40°, we obtain z/a ^ 3.3, fitted to the 15 GHz data, whereas 
SZ07 found z/a =^ 3 J; thus there is a good agreement. On the other 
hand, we find z/a ^ 5.8 using the full model of equation l |16l l. The 
higher value found by us is due to our inclusion of adiabatic cool- 
ing, which reduces the temperature at large radii and thus increases 
the wind opacity there, requiring in turn the emission place at a 
higher height. 

However, although SZ07 used the same total M as that used 
here, they assumed isotropy, without taking into account focusing 
and the reduction of the wind from the irradiated surface of the 
donor. This resulted in a relatively high density of the wind. As a 
result, the 15 GHz emission of Cyg X-1 was found to be attenuated 
to A ^ 0.2, see fig. 13 in SZ07. As also shown there, the average 
attenuation decreased with decreasing frequency at v S 30 GHz, 
to A 0.5 at 2.25 GHz, as well as also decreased with increasing 
frequency at y ^ 30 GHz. However, as argued in Paper I, this is 
not compati ble with the 2-220 GHz spectrum being rather straight 
with ff = jpender et al]|2000t) . and with no substantial dip at 15 
GHz, see also figs. 10-1 1 in Paper I. This supports the wind den- 
sity in the polar regions being substantially lower than that for an 
isotropic wind, see Section |3] Thus, the free-free optical depth of 
the wind is most likely substantially lower than that modelled by 
SZ07. Furthermore, the model of SZ07 implies a very strong or- 
bital modulation, by - 90 per cent in the 1 50-220 GHz range. This 
range was observed bv lFender et al.l(l2000l) , and no apparent orbital 
modulation was seen. In particular, there was one 146 GHz obser- 
vation at ^/2n = 0.95, at which the measured flux was about equal 
to the average one. 

Our best data are those at 15 GHz, and we fit them with the 
model of equation l |16t . We first give the results neglecting the ex- 
tended character of the radio emission, i.e., for b = 0. We obtain 
/ ^ 0.060 ± 0.005, z/a ^ 0.57 ± 0.01, 0nii„/2;r ^ 0.13 ± 0.01, 
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Figure 3. The observed (i.e., at = 0) orbital modulation (data points with 
errros) at (a) 15 GHz, (b) 8.3 GHz, (c) 2.25 GHz fitted by the model (solid 
curve) of equation l |16t with (' = 31°. The flux averages are (F) ^ 11.6 mJy, 
14.2 mJy, 14.5 mJy, respectively. 
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Figure 4. Variations of the free-free opacity along a photon path for our 
model for 15 GHz (at = 0, / = 0.06, ( = 31°). The solid, dashed and 
dotted curves correspond to the orbital phase of {<f) - 0n,in)/2;r = 0, 0.5 and 
1 , respectively. 



at xl - 3.4/8. Fig.|3la) presents a fit of this model to the folded 
and averaged light curve. The model modulation depth is D ^ 
0.30 ± 0.01, and the orbit-averaged attenuation is A = 0.76 + 0.01. 
The values of ^min/ZTr and D are very similar to those obtained by 
iLachowicz et alj ( l2006h . who used a smaller data set and a different 
model. 
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Figure 5. An example of the radial temperature profile for our model for 15 
GHz (at 6 = 0, / = 0.06, ( = 31°) at the direction of & = 45°, = 90°. 
The dashed and dotted curves correspond to the solution of the dift'erential 
energy equation jlOt with the boundary condition imposed at i.2R, and 
1.3R,, respectively, and the solid curve corresponds to the algebraic solution 
of equation U2) . 
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Figure 6. The observed (b = 0) orbital modulation at 15 GHz fitted by the 
model of equation (16) with ! = 31° for the fluxes (a) < (F)/1.15 and (b) 
> 1.15(f), where (F) = 1 1.6 mJy is the average for the entire sample. 
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Figure 7. The dependence of the fitted height along the jet, z (in units of a) 
on the 15 GHz flux, for the case b = 0.5, / = 0. 1. The line shows the best 
fit. 
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Fig. m shows variation of the absorption coefficient along the 
photon path, for our best-fit model. We see the opacity drops fast 
at Z > a. The region / < a is at radii of r ~ 1.5a. Fig. |5] shows an 
example of the radial temperature profile for the best-fit model. We 
see that indeed the solution is virtually independent of the adopted 
inner boundary condition and the differences between the solution 
of the differential energy equation and the approximate one with 
the advection term are relatively small, S 20 per cent. 

We then consider the effect of taking into account the effect of 
a part of the radio flux being emitted well beyond the orbit, assum- 
ing h = 0.5. We obtain the same values for / and 0n,in as for b = 0, 
and z/a - 0.23 ± 0.02 atxl ^ 6.0/8, and with D ^ 0.31 ± 0.01, 
A = 0.82 ± 0.01. We note, however, that the obtained z is rela- 
tively close to the binary plane, which may no more be in the polar 
region, and thus not to be consistent with the fitted low value of 
wind-density scaling factor, /. Thus, we also consider models with 
fixed larger values off, but still « 1. At / = 0.1, z/a 0.62 ±0.02, 
D ^ 0.30 ± 0.01, A = 0.70 ± O.Ol,xl - 13.7/9, and at / = 0.2, 
z/a ^ 1.21 ± 0.02, D ^ 0.29 ± 0.01, A = 0.52 ± Om,xl ^ 27.5/9. 
Since z ~ a in both cases, the emission region is likely to be in 
the low-density polar region. Although these values ofxl are rather 
high, the bad fit quality may be due to our lack of knowledge of the 
actual structure of the stellar wind and of the spatial distribution of 
the radio emission along the jet. Hereafter, we adopt / = 0.1, where 
both / <K 1 and z ~ a, but also use / = 0.2 for comparison in some 
cases. 

Based on general f eatures of the jet emission, e.g., BK79, 
iHeinz & SunvaevI ( l2003ll . we expect that the height at which the 
bulk of emission at a given frequency emerges increases with in- 
creasing radio flux. This is due to the associated increase of the den- 
sity of the relativistic electrons in the jet, which increases the syn- 
chrotron self-absorption optical depth. Assuming the magnetic field 
strength to be proportional to the equipartition valu e with the rela- 
tivisti c electrons, we derive (using the formalism of lZdziarski et al.l 
|20l3). 



ZocF" 



s + 6 
2s + 13 ' 



(32) 



where s is the index of the electron distribution (oc where y 
is the Lorentz factor), which gives q =^ 0.47 for 0.7 < i- < 3.5. 
[Equation (5) of Heinz ^2006 ) gives this result for i- = 2.] The 
increased height of the point where the bulk of emission occurs 
results in less wind absorption and a lower modulation depth. Also, 
a higher X-ray flux (which is positively correlated with the radio 
flux, e.g.. Paper I) leads to a higher temperature, equation (I12t . and 
thus less wind absorption of radio photons. 

Given the good statistical quality of the 15 GHz light curve, 
we can test the above theoretical prediction with the data. For this 
purpose, we have divided the observed range of the 15 GHz flux 
in the hard state, F, into three parts, satisfying (1) F < (F)/1.15, 
(3) F > 1.15(F), and (2) that in between. As the criterion, we have 
used the radio fluxes averaged over a single orbital period, F, in- 
stead of individual measurements, F, since the latter are themselves 
affected by the orbital modulation. The unweighted geometric av- 
erage fluxes in these ranges, j = 1, 2 and 3 are Fj = 8.7 mjy, 
1 1.6 mJy and 15.6 mJy, respectively. In order to account for a part 
of the radio flux emitted at large distances, equation ( 128) , we sub- 
tract from each of the three folded light curves a part of the local 
average, i.e., bFj. 

We take into account that the plasma temperature of the ab- 
sorbing wind changes in response to the changing flux within the 
hard state. Given that the X-ray flux in the hard state of Cyg X-1 is 



correlated with the 15 GHz flux as Fx oc F*''-^, we adjust accord- 
ingly the normalization of Lx and Lion- 

Then, fitting the model of equation l |18t for b = 0, i = 31°, 
/ = 0.060, we have obtained the location of the emission at zj/a = 
0.54 ± 0.03, 0.64 ±0.03, 0.72 ±0.05, which imply = 0.35 ±0.02, 
0.26 ± 0.02 and 0.20 ± 0.02, and fj = 0.32 ± 0.01, 0.26 ± 0.01, 
0.20±0.01, for j = 1, 2, 3, respectively. The folded light curves and 
the fits are shown in Fig.|6] Alb = 0.5, / = 0. 1, Zj/a = 0.59 ± 0.02, 
0.70 ± 0.02, 0.79 ± 0.04, which imply Dj = 0.35 ±0.01, 0.26 ± 0.01 
and 0.20 ± 0.01, and fj = 0.72 ± 0.01, 0.57 ± 0.01, 0.45 ± 0.01, 
for j = 1, 2, 3, respectively. This z(F) dependence is shown in 
Fig. [7] We fit it as a power law, z «- F", and obtain v 0.51 ± 
0.10, i.e., the statistical significance of an increase of z with F is 
5(7. We see that the orbital modulation depth strongly decreases 
with increasing F, also at a high statistical significance. Thus, these 
calculations confirm the dependence of the radio emission region 
along the jet on the jet power in relativistic electrons with a high 
statistical significance. 

As discussed above, we assume that the long term variability 
occurs in the same way for the modulated and unmodulated parts of 
the radio flux. We can then write the dependence of equation ( 132b 
for the intrinsic flux, equation ( 129) . and estimate. 



Z(^iiitr) - ZO 



where zo = Z2, Fo = ^2 + (1 - b) exp(T2)], and 



(33) 



ln(z3/zi) 



ln(F3/Fi)-l-ln 



b + (l - b)exp(T3) ' 
b + il - b)exp(fi) 



(34) 



We obtain q =^ 0.62 ± 0.21, 0.70 ± 0.14 for (b = 0, f = 0.06), 
(b = 0.5, / = 0.1), respectively. In the latter case, we also fit the 
z(Fjn,r) dependence, and obtain q ^ 0.72 ± 0.14, i.e., an almost 
identical value to that of equation ([34}. These values are somewhat 
higher than the theoretical prediction of equation (I32K This may 
reflect either the approximate character of the assumed disc model 
of BK79, the fitted model being too approximate, or both. 

We then fit the 8.3 GHz and 2.25 GHz data, see Fig.[lb-c). 
The obtained parameters at (fe = 0, / = 0.06) are z/a ^ 1.79 ±0.22, 
6.17 ± 1.60 (xl ^ 9.5/9, 3.9/9), which resuU in £) ^ 0.11 ± 0.03, 
0.04 ± 0.03, and A = 0.77 ± 0.04, 0.69 ± 0.13, respectively. The 
phase lags are (p^i„/2n =^ 0.15 ± 0.05, 0.22 ± 0.13, respectively. At 
{b = 0.5, / = 0.1), z/a ^ 1.87 ± 0.16, 6.32 ± 1.16 ^ 18.7/9, 
7.6/9), which result in D ^ 0.11 ± 0.03, 0.04 ± 0.03, and A = 
0.73±0.04,0.64±0.12, respectively. (At / = 0.2, z/a ^ 2.78±0.21, 
8.84± 1.58, respectively.) The phase lags are unchanged. These data 
cover a much shorter interval in the hard state than the Ryle/AMI 
15 GHz data do, and have been taken with a substantially lower 
sensitivity. Thus, the resulting folded light curves are much more 
noisy, especially that for 2.25 GHz, which has the parameters only 
weakly constrained. Fitting in the log space at 2.25 GHz, 8.3 GHz 
and 15 GHz withz cc y-", wefindn 1.4±0.1,;i 1.2±0.1 for/ = 
0.1, 0.2, respectively. This is somewhat more than the prediction of 
the model of BK79, which may be due to either the departures of 
the jet in Cyg X-1 from that model and/or our fitting model being 
too approximate. 
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Figure 8. Comparison of the radio/X-ray correlation f or Cyg X-1 in the h ard state corrected for free-f ree absorpti on (assuming b = 0.5, f = OA, d = 2 kpc; 
black crosses) wi th those for the LMXBs GX 339^ jCorbel et al]|2003l ; green crosses), V404 Cyg l lCorbel et al. 2008; magenta squares) and H1743-322 
iCoriat et al.ll201 It red circles). The dotted and dashed lines have p = 0.6 and p = 1.4, respectively, approximately fitting the two apparant branches of the 
correlation for the three sources. See Paper I for the un-corrected Cyg X-1 data. The Cyg X-1 radio data are for (a) 15 GHz, (b) 8.3 GHz, and (c) 2.25 GHz. 



6 EFFECT OF FREE-FREE ABSORPTION ON THE 
RADIO/X-RAY CORRELATION 

Cyg X-1 in all spectral states shows a power-law coiTelation be- 
tween the radio and X-ray fluxes, F oc F^, where Fx is either the 
flux in an X-ray energy interval, the bolometric flux, Fboi, or the flux 
in the Comptonization X-ray component, see Paper I. In particular, 
the correlation of the 15 GHz flux with Fboi in the hard state has 
the index of pboi = 1 .68 ± 0. 1 1 for the simultaneous ASM and BAT 
data. Qualitatively, the hard-state correlation in Cyg X-1 is similar 
to that seen in the hard state of lo w-mass X-ray bina ri es (LMXBs) 
containing black-holes, see. e.p..[Gallo et al.l j2003h . ICorbel et al.l 



( l200(]|.l200ll2004l2008l) , [Coriat et al.l ( 1201 1[) . However, the corre- 
lation in Cyg X-1 is signi ficantly steeper than th at in LMXBs, e.g., 
Pboi - 0.8 in GX 339-4 fcdziarski et"aLll2004h . A part or whole 
of this difference is due to the effect of free-free absorption of the 
radio emission in the stellar wind. 



An observed correlation index, p = dlnF/dlnFx, where Fx 
is either the flux in an X-ray band, the bolometric flux or the flux 
emitted by hot electrons (see Paper I), is related to its corresponding 
intrinsic index, pi,,,,. = dlnF|„(|./dlnFx, by 

dlnF,„„- dln[fc + (l -fo)exp(f)] 

Pint, = rp, r = = 1 + — — < 1. (35) 

d In F d In F 

So, the effect of removing free-free absorption of the radio fluxes 
is to reduce their variability range, without affecting the variability 
range of the X-ray flux. Note that the factor r is the same regardless 
of the choice of the energy range for the X-ray flux. 

The coefficient r can be determined from the results shown in 
Figs.|6}l7] We estimate it as r ^ \ +Mw[b + (\ - i)exp(f)] /AlnF. 
We obtain r ^ 0.79 ± 0.02, 0.71 ± 0.02 for (b = 0, f = 0.06), 
(b = 0.5, / = 0.1), respectively. This implies pi„ti.boi - 1-33 ± 0.10, 
1.19±010, respectively, for the correlation with the bolometric flux. 
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We note that although the wind absoiption is significantly stronger 
for the b = 0.5 case, it acts only on a 1/2 of the radio flux. 

We can also directly calculate the effect of the wind absorption 
on an individual measurement of the 15 GHz flux, F, performed at 
an orbital phase, (p. Assuming / and / and using the dependence of 
z(^inti) of equation J33l > allows us to write Fj,,,!. as, 

Fi„,, = Fjfc + (l-i)exp[T(z,0-0„„„)]) = Fo|^) , (36) 

which we numerically solve for z, which yields Finn . The resulting 
Fi„ij points for b = 0.5, / = 0.1 are shown together with the 3-9 
keV fluxes in Fig.[8ja). Fig.[8ja) also compares the radi o/X-ray cor- 
relati on for Cyg X-1 with that in three LMXBs (from lCoriat et al.l 
I2OI ih . This figure is similar to one shown in Paper I, which, how- 
ever, showed the Cyg X- 1 radio fluxes as observed. We see that the 
Cyg X-1 points have now the sl ope of p - 1.2, rel atively similar to 
that of the LMXB H1743-322 jCoriat et all 201 1). 

The 8.3 GHz and 2.25 GHz data do not allow us to directly 
make similar estimates, given the poor statistics of the folded light 
curves, see Fig. [3lb-c). S ince the radio spectrum is approximately 
flat, F(v) oc v" ( Ipender e t al. 2000), we assume that the A In F[^ti: 
estimated for 15 GHz is accompanied by changes of z/a at the other 
frequencies proportional to those inferred at 15 GHz, i.e., we use 
equation J33I I with the value of q obtained for 15 GHz, but zo and 
f are from the fits to the 2.25 GHz and 8.3 GHz data. Figs.[8tb-c) 
show the resulting Fj,,,!. points for b = 0.5, / = 0. 1 . Given the strong 
noise in those data, it is difficult to accurately determine the slope 
formed by those points, but it appears similar to that found for the 
15 GHz data. 

However, given the observational uncertainty regarding the ac- 
tual wind density in the polar regions and possible systematic un- 
certainties of the fitted model, the case with a higher wind den- 
sity is also possible, as discussed in Section |5] For example, for 
{b = 0.5, / = 0.2), r = 0.39 + 0.06 at 15 GHz, which corresponds 
to Pintr.boi - 0.65 ± 0.07. Thls would make Cyg X-1 similar to ob- 
jects in the main branch of the radio/X-ray correlation, e.g., GX 
339—4. Thus, our results do not allow us to unambiguously deter- 
mine the slope of the intrinsic radio/X-ray correlation in Cyg X- 1 . 
On the other hand, if we assume that Cyg X-1 in not a unique ob- 
ject, the above constrains the wind density in the polar region to 
/ < 0.2. 



7 ORBITAL MODULATION AS A FUNCTION OF 
FREQUENCY 

Our results for the 2.25-15 GHz data allow us to make some pre- 
dictions for the orbital modulation and attenuation by the wind of 
jet emission at other frequencies, especially higher ones. As esti- 
mated in Section [6] z Z v"' ''-v"'^. Given the implied low values 
of z/a for high-frequency radio emission, it may originate within 
the dense, focused, wind, which would result in much higher wind 
density than in the polar regions, and hence much higher depths of 
orbital modulation than those we estimate here. 

Still, including the focused wind would introduce additional 
free parameters to the model, in particular, the photon path going 
through different wind components. Instead, we assume / = 0.2, 
which, as discussed in Sections|5}l6l approximately corresponds to 
the highest density of the polar wind compatible with the obser- 
vational constraints. Our predictions are shown in Fig. |9] It shows 
the modulation depth (plotted as 1 - D) and the average attenuation 
as functions of frequency. We see that the maximum modulation 
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Figure 9. The orbital modulation depth (plotted as the maximum relative 
attenuation, 1 - Z), solid curve) and the average attenuation (dashed curve) 
as a function of frequency predicted by the model fitted to the 15 GHz data 
with b = 0.5, / = 0.2, i = 31°, and z oc v^' ^ based on the 2.25-15 GHz 
data. 

depth occurs for v 30 GHz, and it is relatively small at higher 
frequencies. Also, the attenuation is relatively moderate, A ^ 0.5, 
at all frequencies. 

We see that we predict much lower modulation depths and 
attenuation than SZ07. The cause for these differences is that our 
models yield a much lower wind density in the absorbing region 
than that assumed by SZ 07, which is supp orted by both the cur- 
rent observational results ( Gies et al. I I2OO8I) , as well as the results 
of our study of the radio/X-ray correlation in Cyg X-1. Our re- 
sults may explain the lack of observed modulation at 146 GHz and 
the 2-220 GHz spectrum being approximately a single power law 
(Fender et al. 2000). Still, given the uncertainty related to the effect 
of the focused wind, it would be of high importance to observation- 
ally study the orbital modulation of high-frequency radiation. 



8 DISCUSSION AND CONCLUSIONS 

We have studied free-free absorption in an optically-thin X-ray ir- 
radiated stellar wind. We have taken into account Compton heating 
and cooling, photoionization heating, bremsstrahlung and line cool- 
ing, and advection. We have derived relatively simple formulae for 
the free-free optical depth, in Section |2] We have also considered 
constraints from eclipses, counter-jets (Appendix |A) and Compton 
scattering (Appendix iBt. 

We have then considered the current constraints on the param- 
eters of t he Cyg X-1 bina ry and its stellar wind (Section |3). As 
found bv lGies etal] ilOOA the wind is strongly anisotropic, with 
the density close to the orbital plane and in the region shadowed 
by the donor much higher than that in the polar regions, which are 
responsible for absorption of the radio emission of the jet. 

In Section O we have discussed the radio data, including the 
available information about the extended character of the emission. 
It appears that a relatively high fraction of the radio emission is 
emitted far away from the binary orbit, and thus it is not orbitally 
modulated by the wind. This effect needs to be taken into account 
in modelling of the modulation. 

In Section [5] we have applied the formalism of Section |2] to 
the available radio data from monitoring of Cyg X-1, taking into 
account the extended character of the emission. In the hard spectral 
state, the radio light curves show significant modulation at the or- 
bital period, which is due free-free absorption in the wind. We have 
fitted the folded and averaged light curves at 15 GHz, 8.3 GHz and 
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2.25 GHz, and obtained the location of the bulk of the emission re- 
gions along the jet. In the units of the separation, our best models 
gavez/a - 0.6-1.2, 1.8-2.8, 6.2-8.8, respectively. The height at 15 
GHz corresponds to ~ 3 x 10'- cm ^ 1.5 x lO^^g, and that at 2.25 
GHz is about 10 times higher. These locations, corresponding to the 
heights at which the jet becomes optically thin to free-free absorp- 
tion, appear to follow a law relatively close to z oc y"', expected 

theoretically (BK79). 

On the other hand. Hein3 ( 1200(f ) in his estimates of the power 
of the jet in Cyg X-1 based on the model of BK79 used the VLBA 
images of SOI, which imply that 50 per cent of the 8.4 GHz emis- 
sion comes from distances ^ 1.5 x lO'"* cm, which is ~ 50a. If 
the entire jet would follow the emission profile of BK79, only a 
tiny fraction of t he emission would take place at z ~ a (see fig. 
I of lHeinzllioO^ . and the resulting orbital modulation would be 
tiny and probably unobservable, whereas it is at the level ~10 per 
cent. Given that the distance of the emission should scale as v"' 
(BK79), it would have been even more difficult to explain the ~30 
per cent modulation of the 15 GHz emission. A likely solution to 
this problem is that the observed resolved emission comes from a 
secondary dissipation event in the jet whereas the emission of the 
inner jet, at z ~ a, is comparable to the resolved one, and it may fol- 
low the model of BK79 . This would a lso solve the problem of that 
jet power calculated bv lHeinzlJioO^ being much lower than that 
inferred from the o bserved large-sca l e optical nebula, pr esumably 
powered by the jet jGallo et al j2005l ; lRussell et al.ll2007l) . Further- 
more, this would solve the problem of the jet base inferred from the 
VLBA measurements being an order of m agnitude higher than that 



theoretically e xpected jRahoui et 



20TT1) . 



Recently, IZdziarski etal] ( l20I2h have calculated the location 



the jet base in Cyg X-1 (defined by the lowest distance at which 
significant emission takes place) at ~ 2 x lO^R,,. U sing the mea- 
surem ent of the turnover frequency of ^ O.I eV of iRahoui et al.l 
( [201 l|), and assuming z oc v ' (BK79), gives a location of the 15 
GHz emission in agreement with our results. 

Our best data, at 15 GHz, allow us to study the orbital mod- 
ulation as a function of the radio flux. We have found a strong de- 
pendence of the modulation depth, changing from D = 0.35 ± 0.01 
to 0.20 ± O.OI when the flux increases by a factor of about 2. This is 
naturally interpreted as the height corresponding to the bulk of the 
emission increasing with the radio flux as a consequence of the cor- 
responding increase of the jet power. In our case, we find z approx- 
imately proportional to the power of about 0.7 ±0.1 of the intrinsic 
radio flux (i.e., that before free-free absorption), which dependence 
is somewhat stronger than the relation in the BK79 model. 

Based on the above results, we have estimated the effect of 
free-free absorption on the form of the hard-state radio/X-ray cor- 
relation in Cyg X-I, which correlation was the subject of Paper I. 
We find the absorption has a substantial effect, but, given the un- 
certainty about the wind density in the polar regions of Cyg X-1, 
we cannot accurately quantify this effect. This effect changes the 
correlation index, p (defined as F oc F^), from about 1.7, found 
in Paper I, to about 0.7-1.3. Thus, Cyg X-1 may have a similar 
correlation to that of the LMXB HI 743-322, which has p ~ lA 
or a number of other LMXBs, e.g. , 



Coriat et al.l 


201lh 


Corbel et al 


2OO3I1 



, GX 339-4 



correlation index is p ~ 0.6-0.8. 

We also discuss predictions for the orbital modulation and at- 
tenuation by free-free absorption of radio emission at frequencies 
higher than 15 GHz. Given the relatively low density of the po- 
lar wind in Cyg X-I, these effects are likely to be much weaker 
than estimated before (in SZ07), which may explain the apparent 



lack of orbital modulation at v = 146 GHz, and the 2-220 GHz ra- 
dio spectrum of Cy g X-I being represented by a single power law 
dFender et al]|2000l) . Still, it would be highly desirable to actually 
measure the modulation depth at high frequency, which would give 
us significant constraints on the structure of both the wind and the 
jet in Cyg X-1. In particular, this would test the possible location of 
the high-frequency radio emission within the focused wind compo- 
nent. 

Finally, we mention that another high-mass X-ray binary with 
strong wind, Cyg X-3 (which is likely to contain a black-hole) does 
not show any orbital modulation of its very strong radio emission. 
As discussed by SZ07, this may be due to the binary orbit of that 
system being ~I0 times smaller than that in Cyg X-I. If the jet of 
Cyg X-3 is similar in size to that of Cyg X-I, this would imply 
values of z/a ~10 times higher, which is likely to reduce the depth 
of the radio modulation to unobservable values. 
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APPENDIX A: ECLIPSES AND THE COUNTER-JET 

The jet emission can be obscured by the donor provided z < R^- 
The lack of such eclipses implies, 



; < 7t/2 - arctan — h arccos 



(a2+z2)l/2' 



z<R.. 



(Al) 



In case of the lack of X-ray eclipses (z = 0), i < arccos(i?,/a), 
which implies no eclipses of the jet emission. Thus, equation jAlb 
is of importance only in systems w here X-ray eclipses are observed 
(e.g., M33 X-7. IOroszet ai]|2007l) . in which case observations of 



eclipses of the radio emission or their lack would impose a con- 
straint on the location of the radio-emitting region. In the case of 
Cyg X-1, the lack of X-ray eclipses implies i < 64° for a/R, ^ 2.3. 

There will be no eclipse of the counter-jet emission by the 
donor for either. 



a R, 

i > arctan h arcsm — tttt, 

kl (a^+zy- 



\z\ > R„ 



I < arctan ■ 



R, 



(A2) 



(A3) 



kl {a^ + z^y^ ' 

For example, for a = 2.3R, and z = -a, no counter-jet eclipses 
correspond to ; < 28° or ; > 63°. In addition, an accretion disc with 
the outer radius, r„, will block the jet emission for i < arctan(ro/z). 

Equation ( |18l > also describes the attenuation of the counter-jet 
emission (in the constant velocity, constant temperature approxi- 
mation), with the substitution of z — > -z- Note that then there is a 
singularity in equation ( 118) is at = and z/a = - cot;, which is 
a true singularity, caused by photon paths going through the centre 
of the donor. 

Note that the counter-jet emission is both relativistically 
beamed away from the line of sight and more absorbed than the 
direct jet emission. Also, the counter-jet emission can be obscured 
by the star at some combinations of z and see equation iA2i . 
All these factors can account for the lack of detection of the radio 
counter-jet in Cyg X-1 (SOI). 



APPENDIX B: COMPTON SCATTERING AND 
BOUND-FREE ABSORPTION 

We now consider scattering or absorption proportional to the den- 
sity (which is the case for bound-free absorption) of photons emit- 
ted by the jet and accretion flow. In the latter case, we require 
that the medium temperature and ionization parameter are constant 
along the photon length, which is a commonly adopted assumption 
in calculating bound-free absorption. 
The optical depth is then. 



(Bl) 



where, in the case of Compton scattering, Tq = a-cnoa/i Jii^, and 
(Tc is the scattering cross section. We can integrate equation ( IBll ) 
numerically, or, assuming v = Voo, we obtain. 



r n/2 + arctan(i-M ''^) 
To u 



1/2 



(B2) 



As before, m = for (p = k, z/a = cot /, and equation jB2b appears 
to have a singularity. However, r is finite, t/tq = sin The other 
singularity at = 0, z/a = - cot /, is a true one. At z = 0, equation 
( [Bit yields 

At 2i 

— = :. (B3) 

To cos I 

Cyg X-1 shows an orbital modulation at hard X-rays, at which 
energies the bound-free absorption is negligible, and which is most 
likely due to Compton scattering of photons away from the line of 
sight by the wind electrons. The source of X-rays in this case is 
close to the compact object, and thus we may set z = 0. However, 
the assumption of v = Voo leads to a significant underestimate of 
the actual t. Thus, we use the numerical integration of equation 
dBlt . For the parameters of Cyg X-1 listed above, we obtain To - 
8.4 X 10 ~^/, where the Thomson cro ss section has been used . The 
results of iLachowicz et all ( l2006h and lZdziarski et all ( l2011al) yield 
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the orbital modulation depth at ~20-100 keV hard X-rays (which 
photon interactions are dominated by the Thomson scattering) of 
0.03-0.04. We find D ^ 0.02/ at i = 31°, which implies / ^ 1.5- 
2. This appears to account for the wind density being enhanced near 
the orbital plane (see Section[3]l. 



